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STABILITY DERIVATIVES OF CONES AT SUPERSONIC SPEEDS 
By Murray Tobak and ~Wi1.I1.am R. Webrend 


SUMMARY 


The aerodynamic stability derivatives due to pitching velocity and 
vertical acceleration are calculated by use of potential theory for cir- 
cular cones traveling at supersonic speeds. The analysis is based on two 
theoretical techniques used successfully previously in application to the 
case of uniform axial and inclined flow. In the first, potential solu- 
tions for axial flow and crossflow are derived from the first-order wave 
equation but in application to calculations for the forces no approxima- 
tions are made either to the tangency condition or to the isentropic 
pressure relation. The second method consists in combining the first- 
order crossflow potential with an axial- flow potential correct to second 
order. Closed-form solutions by both methods are found for a cone,' and 
numerical results for the stability derivatives are presented as a func- 
tion of Mach number for cones having semivertex angles of 10° and 20°. 

In addition, expressions for the forces, moments, and stability 
derivatives of arbitrary bodies of revolution are obtained using Newtonian 
impact theory. Numerical results for cones compare well with those 
obtained from the combined first- and second-order potential theory at 
the highest Mach number for which the latter theory is applicable. 


INTRODUCTION 


The importance of the body as a lift -producing component of aircraft 
flying at supersonic speed has occasioned a great deal of theoretical work 
from which it is now possible to calculate the body* s static aerodynamic 
properties with good accuracy (see, e.g., refs. 1 to 3 and attendant 
bibliographies). There is required, however, along with the static prop- 
erties, theoretical information from which the dynamic behavior of bodies 
can be calculated, and in this field no work has been done that can be 
said to be applicable to nonslender bodies traveling at high supersonic 
Mach numbers. Slender-body theory, as is well known, fails to predict a 
dependence of the aerodynamic coefficients either on Mach number or on 
body shape (see, e.g., ref. 4). The work of Dorrance (ref. 5 ) based on 
the linear theory does indicate a Mach number dependence, but here too 
approximations made in the analysis effectively limit its application to 
bodies of vanishingly small thickness. 
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In this report, an attempt is made to overcome these limit ations "by 
adapting to the calculation of the body*s rotary stability derivatives two 
theoretical methods derived by Van Dyke (ref. l) that have proven suc- 
cessful in calculations for the static aerodynamic derivatives. It was 
shown by Van Dyke that the solution for normal force of a cone derived 
from the first-order potential equation may be greatly improved in accu- 
racy (in comparison with the exact numerical results, ref. 2) if no 
approximations are made to the tangency condition or to the isentrcrpic 
pressure relation. It was also shown that a further improvement could 
be realized by the use of a combination of first- and second-order poten- 
tial solutions. The same ideas are used herein to calculate for a cone 
the stability derivatives due to pitching velocity, and C^, and due 

to vertical acceleration, Cj^. and C^.. The results are believed to be of 

the same order of accuracy as those of reference 1 for the normal force . 
due to angle of attack, and hence, for a cone of given thickness ratio, 
should apply to the same range of Mach numbers over which the normal force 
compares well with the Kopal results (ref. 2 ). Moreover, as in the angle- 
of -attack case, the cone solutions given herein are adaptable to the cal- 
culation of the stability derivatives of other more general body shapes 
by use of the techniques described in references 6 and J. 


In addition, expressions for the forces, moments, and stability 
derivatives of arbitrary bodies of revolution are obtained from Newtonian 
impact theory in order to furnish some information about the nature of 
these quantities at Mach numbers beyond the highest for which the results 
derived from potential theory are applicable. 


NOTATION 


a 0 

°N 

c m 


speed of sound in still air 

normal force 


normal-force coefficient, 


<3o S 


pitching-moment coefficient, ffcraent 

9o SZ 


pressure coefficient 


axial-force coefficient. 


axial force 


*o S 


body length 


Mach number, ^ 
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q angular velocity (sketch (h)) 

q Q dynamic pressure, i p Q V 2 

R(x) radius of "body in r,u plane (sketch (a)) 

5 body base area, j® 2 (Z) 

t time 

u,v,w free-stream velocity components relative to cylindrical coordinates 
fixed in body 

x,r,u cylindrical coordinates (sketch (a)) 

V velocity in axial direction of body 

V-jj component of free-stream velocity normal to body surface 
a angle of attack 

0 n/M 2 - 1 

y ratio of specific heats 

6 slope of body meridian curve (sketch (a)) 

p Q density of still air 

t slope of cone surface, 

$ total potential 

<p perturbation potential 

q> Q first-order uniform axial flow potential 
X Q second-order uniform axial flow potential 

ft free-stream potential 

When a, d, and q are used as subscripts, a dimensionless derivative 
is indicated, and this derivative is evaluated as the independent variable 
(a, d, or q) approaches zero and all other variables are identically zero. 
Thus, 
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a=q=o 


a=q =0 


w 


a.-a=o 


AMETSXS 


Coordinate System and Definition of Free Stream 


In the succeeding analysis we consider a pointed body of revolution 
flying at constant supersonic forward speed. Our purpose is to calculate 
the body* s aerodynamic stability derivatives corresponding to the follow- 
ing motions: (l) sinking with uniform vertical velocity, ( 2 ) flying in 

a circular path with uni form angular velocity and at zero angle of attack, 
(3) sinking with uni form vertical acceleration. In order to define the 

A r ,«=7X/2 motions conveniently, a cylindri- 

e cal coordinate system is chOBen 

✓ >s / r - 1 that is fixed with respect to the 

[ ^ body. As shown in sketch (a) , 

( /ii \ x the origin of the coordinate 

T + t ) W:i o N. I * system is placed at the body 

V J r=R(x) apex. The positive branch of 

" J — — J the x axis is coincident with 

the body’s axis of revolution, 
and the coordinates r and 10 

Sketch (a) are measured in a plane perpen- 

dicular to the x axis. With 

respect to this system of coordinates, the components of free- stream 
velocity u,v,w for the three motions are given below. 

Sinking with uniform vertical velocity .- This case is of course 
equivalent to that of a stationary body situated in a uniform inclined 
stream and has already been treated extensively (refs. 1 to 3)* is 
included here again since the methods to be discussed subsequently for 
calculating the stability derivatives due to pitching velocity and verti- 
cal acceleration are in large part derived from the one used here. Let 
a be the angle of inclination of the stream with respect to the to = 0 
plane; then the components of stream velocity in the axial, radial, and 
azimuthal directions, respectively, may be written , 1 


1 Note that we have chosen to designate as V the axial component 
of velocity rather than the resultant flow velocity. It is necessary 
to do this in the subsequent cases, and we comply here for the sake of 
consistency. Note also that the radial and azimuthal velocities are c, 
measured in directions respectively normal and tangential to the body c - 
surface in the r,w plane. 
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-u = V 

v = V tan. a sin w 
w = V tan a cos w 


( 1 ) 


Once the normal-force coefficient Cjj and pitching-moment coefficient 
C m corresponding to this motion have "been calculated, the stability 
derivatives and are formed according to the definitions,- 






a. 


Pitching with uniform angular velocity .- For this case, the body is 
considered to be flying in a circular path at zero angle of attack and 
•with constant angular velocity q. The motion is, of course, that of the 
whirling-arm experiment, and, in terms of the latter case, we specify 
that the point of attachment of the body to the arm be at the body nose. 
With respect to a fixed system of coordinates, the pertinent physical 
dimensions are as shown in sketch (b) . 

It will be assumed that the pitching 
rate q is small, so that ql « V. 

The radius of the flight path V/q, 
is then large compared with the body 
length. The body is assumed to have 
traveled far enough to have outrun 
its starting sound waves, but not so 
far as to have encountered its own 
wake, so that the flow may be said 
to be steady. 

In the body system of coordi- 
nates, the components of stream 
velocity are. 


u = V - qr sin u 
v » qx sin oi 
0 w =* qx COS U) 


( 2 ) 



The stability derivatives due to 
the motion are calculated according 
to the following definitions: 


Sketch (b) 
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Uniform vertical acceleration . 
Since the body sinks (-without 
pitching) at a linearly increasing 
rate, the flight path describes the 
parabolic arc shown in sketch (c) . 
It is assumed that the motion has 
continued long enough so that tran- 
sient effects have died out and, 
further, we choose to begin record- 
ing time at the point where the 
angle of attack is zero. 

As viewed from the body- 
coordinate system, the components 
of free-stream velocity are, 


u = V 

v = Vdt sin u 
w = Vdt cos a) 


( 3 ) 


where d is a constant. Note that although it has been assumed that 
transient effects associated with the start of the motion have disappeared 
nevertheless, the flow in this case is unsteady, the crossflow velocities 
being linear functions of time. 

The stability derivatives due to the motion are evaluated when the 
angle of attack is zero (t = 0) according to: 



Equations of Motion and Boundary Conditions 


For the main part of this paper, the analysis will be based on the 
first-order steady and time-dependent potential equations for compressible 
flow. In body coordinates, the latter equation is. 
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o2* * *r, 2 M, . 1 _ _ n 

S’rr "72 F + a^ ^xt + ^-g *tt _ 0 


A perturbation potential cp is introduced such, that 

$(x,r,w,t) »Q(x,r,w,t) + cp(x,r,u,t) (5) 

f 

where fl is the free- stream potential. Direct substitution gives for J 
the first-order time-dependent perturbation equation, 

< Prb + ^'Ptt =0 ( 6 > 

The steady-state counterpart of equation (6) is obtained by eliminating 
time derivatives, so that, 

( 7 ) 

Equation (7) will be used to determine solutions for the stability deriva- 
r tives due to steady sinking and steady pitching, whereas the stability 
derivatives due to vertical acceleration are derived from equation (6) . 

In either case, however, the boundary conditions have the same form, 
namely, that, 

(l) Velocity perturbations vanish on the Mach cone emanating from 
the body nose: 


<p(x, f , «,t) = 0 


(2) The flow velocity normal to the body surface is identically zero: 


SR _ + flr(x,R,h),t) ^ 

Ax cp x (x,R,u>,t) + fi x (x,R,u,t) 

The velocities % and fix are, of course, the free- stream radial and 
axial velocity components, respectively, so that more conveniently. 


an <P r (x,R,w,t) + v(x,R,u»,t) 

— S 1 — .. 

dx cp x (x,R,u),t) + u(x,R,co,t) 


( 10 ) 
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where, for the three motions to he considered, u and v are given by 
equations (l) to ( 3 ) . Note also that we designate the independent radial 
coordinate by r, whereas in evaluation of conditions at the body surface, 
the fact that r is a function of x is indicated by the use of E (i.e., 
r = R = R(x) ) . 


First-Order Solutions for the Potential 


We consider next the task of finding solutions for the perturbation 
potential satisfying equations (6) or ( 7 ) and compatible with the boundary 
conditions, equations (8) and (9), corresponding to each of the three 
motions . 

Sinking with uniform vertical velocity .- As mentioned previously, 
the method of solution in this case is well known j however, since the 
methods to be used for the two subsequent cases derive from the one used 
here, a brief account of the essential steps is given below. 

The steady perturbation potential is broken into two parts: A poten- 

tial cp Q independent of u and hence corresponding to a uniform axial 
flow , and a crossflow potential cp 1 . The total perturbation potential is 
then the sum of (p Q and cp^ For the axial-flow potential cp Q , the equa- 
tion of motion is, 

'Po rr + ^-»% ac =° 

with the boundary conditions (from eqs. (l), (8), and (lO)) 



cPo^K) 

— ss ... 

dx qp 0x (x,R) + V 


(13) 


A solution to equation (ll) that automatically satisfies (12) is (ref. 6), 




f(x- pr cosh u)du 


cosh 


-l 


x 

pr 


m 
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where f(x- pr cosh u) represents the distribution of sources along the 
x axis, and is to be chosen such that equation (l^) satisfies the 
boundary condition (13) . 

For the crossflow potential cp 1 , the equation of motion is, 

^rr + ^ + ^” P ^xx =0 to) 

with the boundary conditions (eqs. (l), (8), and (10)), 



(16) 


dR _ 
dx 


<P lr U,R,w) 


+ V tan a. sin w 


cp lx (x,R,u) 


(17) 


A solution to equation (15) that automatically satisfies (l6) is (ref. 6), 


o 

(p 1 (x,r,u) = -p sin u J m(x- pr cosh u)cosh u du (18) 

cosh” 1 iL 
pr 

where m(x- pr cosh u) represents the distribution of doublets along the 
x axis, and is to be chosen such that equation (18) satisfies the boundary 
condition (17) . 

Pitching with uniform angular velocity .- The procedure in this case 
parallels the above development. Again, the perturbation potential cp 
is broken into a uniform axial flow potential cp Q and a potential cp 2 , 
each having the same equation of motion as equations (ll) and (15) , 
respectively. Note, however, that Tin like the previous example, the cal- 
culation for the potential cp 2 must take into account the nonuniform 
axial component of stream velocity, -qr sin u> (eq. (2)). The boundary 
conditions on <p Q are, from equations (2), (8), and (lO), 


'Po(*-f)“° to) 

gR_ V X > R) 

<p 0x to R ) + v 


(20) 
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while the boundary conditions on cp 2 are. 



( 21 ) 


,tr + qx sin to 

S ----- (22) 

cp 2x (x,R,o») - qR sin w 

The solutions for the potentials cp Q and cp 2 of course have the same form 
as equations (l4) and (l8) , respectively, where f(x - pr cosh u) and 
m(x-pr cosh u) are to be chosen such that equations (20) and. (22) are 
satisfied. 

Sinking with uniform vertical acceleration .- As mentioned previously, 
this is a problem in unsteady flow, since the crossflow velocities v 
and w are functions of time. It is still possible, however, to consider 
the potential <p in two parts; an axial component cp 0 that is independent 
of both t and w, and. hence is again governed by equation (ll), and a 
crossflow component q> 3 , governed in this case by equation (6) , with 
9 = C P3* 


For the axial-flow potential, we proceed exactly as in the two 
previous cases. The boundary conditions to be satisfied are, from 
equations (3), (8), and (10), 



(23) 


aB _ yUPO 

ax v+ cp 0x (x,R) 


(24) 


The solution for tp o is given by equation (l4) where f (x - (3r cosh u) 
is chosen such that equation (24) is satisfied. 

For the -unsteady crossflow potential qp 3 , we adapt a concept used 
previously in wing theory by Ribner and Malvestuto (ref. 8) and origi- 
nated by C. S. Gardner. It Is easily verified by substitution into 
equation (6) that the following relation satisfies the unsteady potential 
equation no matter what the constant K may be: 

*Tp = Eijr(x,r ,w) + ^t - ^^X(x,r,u) 


(25) 



MCA TN 3788 


11 


where 

tJt steady- state potential for unit pitching velocity about body apex, 

<1 

m 

X steady-state potential for unit angle of attack, 1 — 

Note that the potentials involved in equation ( 25 ) are steady-state 
potentials and, further, are just the potentials cp x and <p 2 developed 
in the two previous sections. Having found a solution for the potential, 
we then must satisfy the boundary conditions which, from equations ( 3 ) , 
(8) , and (10) , are 



( 26 ) 


dR 

dx 


<p^(x,R,w,t) + Vat sin w 
cp 3x (x,R,a>,t) 


(27) 


Equation (26) is satisfied immediately, since the potentials i|r and X 
individually satisfy it. The tangeney condition, equation ( 27 ), is then 
readily satisfied by proper choice of the constant K. 


Pressure, Force, and Moment Coefficients 


It will have been noticed in the previous section that we have speci- 
fied the exact form of the tangeney condition for each of the motions 
considered. The same will be done for the pressure relations, to be given 
below. Following Van Dyke (ref. l) , the view taken on this point is s imply 
that approximations to the tangeney and pressure relations, while justifi- 
able mathematically on an order baBis, serve to impair unnecessarily the 
accuracy of the solution in comparison with known numerical results. While 
this is known to be the case only for the -uniformly sinking motion, the 
formulation of the problem for the other motions has introduced no further 
approximations, and hence it is entirely reasonable to suppose that the 
same order of accuracy will be realized in these cases as well if~ the 
pressure and tangeney relations are not approximated. 

The pressure coefficient must be defined separately for each of the 
three motions^ this has been done in the appendix, and the results are 
repeated here for convenience. 



12 


NACA TR 3788 


Sinking with uniform vertical velocity . - 


7 M 2 


7 

, 7=1 


(l- - 1 


(28) 


where. 


A 0 . Aa <Pr \ c Px 2 + c Pr 2 + ( < Po)/ r ) : 

Ai=— + 2 tan o(-^ eos <o+ — sin ^ 


and r is taken as R(x) when evaluating Cp at the body surface. 
Pitching with uniform angular velocity . - 

where 

. ^Px 2cp x 2qx <Pr . \ fox 2 +9r 2 + ( ( Pw/ r ) 2 ~l 

V V 2 V \rV V J L V 2 J 

Sinking with uniform vertical acceleration . - 




7 _ 

7-i 


-1 


(29) 



Wormal-force and pitching-moment coefficients .- Once having the 
pressure coefficient, one can determine the normal-force and pitching- 
moment coefficients from the following relations: 


Cjj = J' R(x)dx J' Cp(x,R,a>)sin u dto 


(31) 
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o 


+ tan 0R e 


(x)Jdx J Cp(x,R,t»))sin w du 


(32) 


where the pitching moment is referred to an axis through the body nose. 


Stability Derivatives for a Cone - First-Order Theory 


The foregoing equations for potential and pressure coefficient axe, 
in principle, applicable to the determination of forces, moments, and. 
stability derivatives of arbitrarily shaped bodies of revolution. In 
order to complete the analysis, it is necessary to find distributions of 
sources and doublets f(x- pr cosh u) and m(x - pr cosh u) that are com- 
patible with the boundary conditions corresponding to the specified body 
shape. Unfortunately, analytical expressions for these quantities have 
been fo und only for the cone. However, it has been shown in the axial- 
and inclined-flow cases (refs. 6 and 7) how, by the use of summation 
techniques, the cone solution can be used to find solutions for other 
more general body shapes to any desired accuracy. The same techniques 
are readily adaptable to the other motions considered herein. Therefore, 
as a necessary beginning toward obtaining the stability derivatives of 
more general body shapes, the calculations for the cone are carried out 
below for each of the three motions under consideration. 

Sinking with uniform vertical velocity .- Consider first the potential 
for axial flow cp Q , given by equation (14) . For the case of a cone, whose 
surface is given by R = rx, an appropriate distribution of sources is 
simply (ref. 6), 

\ 

f (x - pr cosh u) = Aq(x - pr cosh u) (33) 

where A 0 is a constant, to be determined from the boundary condition, 
equation (13) . Integrating equation (l4) , and substituting the appro- 
priate derivatives in equation (13) , we get for A 0 , 


A-o - 


Vt‘ 


T^COBh " 1 — + *Vl- p 2 -! -2 
pT 


(34) 


whereupon the axial- flow potential is. 


<P Q (x,r) =Ad^/ X 2- p 2 r2- x cosh" 1 ^ 


(35) 
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The computation for the crossflow potential <Pi proceeds in the same 
way. An appropriate distribution of doublets is found to be (ref. 6 ) 

m(x - pr cosh u) = Ai(x - pr cosh u) ( 36 ) 

Substituting in equation (l8) , integrating, and solving for Ax in 
equation ( 17 ) , 


A, = 


2VT 2 tan a 

(2-r 2 + 1) Jl - p 2 !* 2 + p 2 r 2 cosh“ 1 — 

Pt 


(37) 


and 


<pi(x,r,w) = Ai ^ sin ui 


(tJ* 


P 2 !* 2 - Pr cosh' 


- 1 — ) 
Pr/ 


( 38 ) 


The perturbation potential cp is then given by the sum of equations (35.) 
and ( 38 ). 


Having obtained the potential, one can now compute the pressure 
coefficient by equation ( 28 ) . However, it is more convenient to proceed 
directly to the calculation of the stability derivatives Cjj a and C^, 
by means of the expressions. 




R(x)dx 



■\ 


> 


Cma 



_2_ 

S2 


^[xR(x) + tan 0 R^xJldx J' 

o Jt a 


sin <0 dm 
o 


( 39 ) 


where 
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("£*) “ - i 1 - *^ (Al) a -» o] 

a o 


7-1/ 


CG 


The quantity 


(' 


d°Pj\ 


da 7 

a — > o 


is computed by substituting the appropriate 


derivatives of <p into equation ( 28 ) and differentiating with respect 
to a as indicated. Substituting the results into equation (39) and 
integrating gives (letting y =' 1.4), 

*\ 

a2t 


where 


_ ( 1+tS } 
N a \A + 0 2 / VL + A + 2 tV 


**-- 1(1 + 


P 2 T 2 cosh " 1 
A = — & T - 

Ji - p 2 t 2 


(40) 


Fl = {l + 0.2^[l - (1 + 


The above results reduce to the slender-body-theory result as the thick- 
ness parameter r approaches zero: 


Cw = 2 
A ’a 

Cm a 



(41) 


* 


Pitching with uniform angular velocity .- Again, consider first the 
uniform axial flow potential cp 0 . There 'is nothing new to compute, 
however, since the equation of motion and boundary conditions are the 
same for all three cases. The potential cp 0 is, therefore, given by 
equation (35) • 
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For the potential cp 2 (eq.. (l8)) a distribution of doublets that 
satisfies the boundary condition, equation (22), is, 

m(x - |3r cosh u) = B^x - f3r cosh u) a (42) 

Substituting equation (42) in (l8), integrating, and solving for Bg in 
equation (22), we get. 


(p 2 (x,r,co) = Bssin u 



where 


1 


cosh -1 

3r 


m 


B2 = 


-qr 2 (l + t 2 ) 

(t 2 - l^^OOSh- 1 i -Jl - P 2 T 2 J - | (1 - 


As before, the total perturbation potential cp is_the sum of tp Q 
(e<l- (35)) and <p 2 (eq[. (43)). Also, in the same way, the stability 
derivatives Cjj and Cm are computed according to. 



where Cp 2 is given by equation (29). There results. 
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C w = — T -fl 1 l ( p2 \ i 1 3,(1- A) - 2(1- 1 1 

% 3 1 l U+P 2 / ^ U(l- ps-r 2 )- 3(1- t2)(1-A)J J 


> (45) 


% - - f a + T 2 )c Nf 


where A and have been defined in equations (40). 


As in the previous case, these results reduce to the slender-body- 
theory result as t -> Oj 



(46) 


Sinking with uniform vertical acceleration .- The axial-flow poten- 
tial cp 0 is again given by equation (35)* For the time-dependent cross- 
flow potential <p 3 , use is made of equation (25). Inserting the values 
for cp! and cp 2 obtained in the two previous cases into equation ( 25 ) 
gives. 



(47) 


where B 2 /q. and A 2 / tan a have been defined in equations (k3) nnfl ( 37 ). 
The constant K is then determined by substituting the appropriate 
derivatives of cp 3 into the boundary condition, equation ( 27 ). There 
re stilts , 
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K = 


M 2 

3 2 (l + T 2 ) 


(3'r 2 + l)Jl^W^ + (1 - T 2 )p 2 T 2 cosh“ 1 ~- 
PT 

(2t 2 + 1)^1 r e. 2 T £ + P^COSh' 1 — 


(W) 


As before, tbe sum of qp 0 (eq.. (35)) and cp 3 (eq.. (47)) is then the total 
perturbation potential 9 due to uniform vertical acceleration. 

The stability derivatives Cw. and Cm. are formed in the same way 

a a 

as the two preceding e xamples , using the pressure coefficient, equa- 
tion ( 30 ) , and the relations. 



I / E(x)ax / 

O £ 

’"2 





\ (49) 




+tan0 


'& ->o 


R 2 (x)]dx / 



sin u du 


d->o 


It should be noted that the stability derivatives are evaluated when the 
angle of attack is zero. This occurs when t = 0, so that, after taking 
the required time derivative in equation ( 30 ) , the remaining terms multi- 
plied by t may be eliminated. The results are. 
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°Hd“l r 4(^){[ 1+ 5— 'M 


(l-A)T 3(1 - a) - 2(1- p 2 -!- 2 ) 


i+a+2t 2 JU(i- p 2 -^) - 3(1 - a)(i- t 2 ) 


]- 


(1 - A) 
1+A + 2t‘ 


1 + JL^L\ 

v i+a+2t y 


>(50) 


°“4 = ' ! U + 


Ab t — > 0, the above results reduce to the slender-body- theory result, 

(51) 


Off- - I 
a 3 


Cm. = - ± 

a 2 


Stability Derivatives for a Cone - Combined. 
First- and Second-Order Theory 


It has been shown by Van Dyke (ref. l) that a further improvement 
in the first-order results for pressure and normal force due to inclined 
flow past a cone can be realized by making a second-order correction to 
the axial-flow potential. This is called in reference 1 a "hybrid theory." 
The same idea can be incorporated in the other two cases considered herein, 
since in all three, the axial-flow potential has the same form. Moreover, 
as in the purely first-order case, the resulting solutions are adaptable 
to the calculation of the stability derivatives of bodies of general shape 
by use of the technique outlined in reference 7» Unlike the inclined-flow 
case, however , it is not assured that by making the second-order correc- 
tion the solutions corresponding to uniform pitching velocity and unif orm 
vertical acceleration are necessarily improved. In the absence of exact 
numerical results with which to compare the solutions, it can only be 
supposed that such an improvement is likely, again in view of the fact 
that in the formulation of the problem no approximations are made beyond 
those which also exist in the inclined-flow solution. There follows a 
brief description of the method of applying the second-order correction 
to the potential for the cone. 

For the case of uniform axial flow past a cone, the potential X Q , 
correct to second order, as obtained from. reference 1 is. 
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Xo 

Vx 



- cosh -1 ~ 
pr 


+ 




E 2 !! 2 


C—*- 1 f) S - tH + «A - W 


cosh” 1 

pr 




J 


where. 


D = 


(H-P-Q)t 2 


J 1 - P 2 i^ + i^cosh -1 “ 



(N-l)p 2 r 2 A _ (H + 1)P 4 t 4 _ 3 P 6 t 4 

(1- ps-r 2 ) (1 - pa-r 2 ) ^ (p 2 + A)(l - p 2 -^) 


f M f _op 2 ft , (W + ljpV 2 | (if - l)p 4 T 2 A [ 1 p 6 (2 + pgr 2 ) 

V^+S 2 / (l-? 2 ! 2 ) ( 1 - fs 2 ! 2 ) 4 ( S 2 +a)(i-^t 2 ) 


E = 


” P 2 t 2 + T-gsech-^-pT 


_ (7 + ipM 2 
2p 2 
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This result applies to the uniform axial-flow potential for each of the 
three motions . The derivation of pressure coefficient and stability 
derivatives then proceeds exactly as before, with the exception that 
equation ( 52 ) is used in place of equation (35) • The results for the 
stability derivatives are. 


°*a = 2 T ^T^) (l AV+d) 




■ I + KaTf ) (1 + 


3(1 -A) - 2(1- 3 2 t 2 ) 


Ml- ^T 2 )- 3(1- a)(i- t 2 ) 


■]} 





t 2 (1-A)T 3(l-A) - 2(1-P 2 t 2 ) 

1+A+2t 2 JU(1- P^t 2 ) — 3 (l**A)(l — t 2 ) _ 



( 1 -A) ] 
1 +A + 2 T 2 J 


(1 - A) 

1 +A + 2 T 2 


St - - f <1 - 

Si - - | (1 + t 2 )C N(1 

St - - i * ^St 


where 


and 
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Note that these results are derivable from the purely first-order solutions 
(eqs . (40), (45), and ( 50 )) if the factor p 2 /(A+p 2 ) is multiplied by n 
wherever it appears in the first-order solutions. Therefore, since |i 
goes to 1 as r 0, these results also reduce to the slender-body-theory 
results (eqs. (4l), (46), and (51)) a- 8 t 0 . Note also in equa- 
tions (53) that the pitching-moment coefficients are simple multiples 
of their respective normal- force coefficients. For a given thickness 
ratio, the centers of pressure are therefore invariant with Mach number. 


Transfer of Axes 


The results for Cjj^, C^, C^, and presented in the preceding 

sections are applicable only to the case of a body whose center of gravity 
is at the nose. These results may be used, however, to calculate the sta- 
bility derivatives for any other cent er-of -gravity position, by means of 
the tr ans fer relations given below. 


Xo 

“ Cm - — C„ 


“'Nr. “ °N 


<3o l 


Xq 


Cm a ~ 0m a o + — C N a 


Sl " So + T %> " ~l ^"“O " \l) 


> 


(54) 


^ - V - t J 

Here, Xo is the new center of gravity position, measured positive rearward 
from the nose, and the subscripted terms axe the stability derivatives as 
calculated for Xq = 0. 


Newtonian Impact Theory 


Kesults canno t be obt aine d by the potential theories used in preceding 
sections above a Mach number for which the Mach cone from the body apex 
coincides with the body surface. Far bodies of moderate thickness ratio, 
this condition limits the applicability of the theories to the Mach number 
range below about 3 or 4. For Mach numbers very much higher than this 
limi t, it may be expected that the underlying assumptions of the Newto ni a n 
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impact theory become increasingly valid. Therefore, in order to provide 
some information on the nature of the aerod yn a mi c forces, moments, and 
stability derivatives at very high Mach numbers, the Newtonian theory is 
adapted below to the derivation of these quantities for arbitrary bodies 
of revolution. 

Pressure, force, and moment coefficients .- The assumption basic to 
impact theory is that the flow, upon striking the body, loses entirely 
its component of momentum normal to the surface and continues along the 
surface with its tangential component of momentum unchanged. The loss 
of momentum normal to the body surface yields a pressure force that is 
simply (ref. 9 ), 

p-p 0 g% a 

^ V 2 

where Vjj is the component of velocity normal to the body. It is to be 
noted that only those portions of the body that "see the flow," that is, 
are under direct impact from the stream, will experience a pressure force. 
The remainder of the body is generally assumed to be at the pressure of 
the free stream, so that on these portions, the pressure coefficient is 
zero. Therefore, when integrating the local pressures over the body sur- 
face to obtain total forces and moments, the integration proceeds only 
over the portion of the body receiving compression flow. 2 The formulas 
for normal force, axial force, and pitching-moment coefficient then are. 



r *u «<%(*) 

C N = “ S J R ( x ) dx J Cpsin w 




dill 


p X u n“u( X ) 

C x = j| / R(x)tan 0 dx J Cpdw 


\ (56) 


Uu(x) 


r> w uv~/ 

C m ~ SlJ t(x " x o) R ( x ) + tan 0 R 2 (x)]dx J Cpsin w dw 


2 Actually, in addition to the impact forces, there should be consid- 
ered the centrifugal forces which arise from the fact that the flow parti- 
cles follow curved paths along the body after impact (ref. 9). Some gross 
estimates of their magnitude have been made which indicate that they axe 
negligibly small for the small changes in a, 4, and q. that are of primary 
interest here. A more precise estimation of their m agnitude is a matter of 
B o r ne difficulty, and in view of the approximate nature of the entire theory, 
the effort required in this direction does not appear to be warranted. The 
centrifugal forces are therefore neglected in the subsequent analysis 0 
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Sketch (e) 


■where the quantities (^(x), 
x u , and Xq are illustrated in 
sketch (d). Provision has been 
made in equation (56) for an 
arbitrary axis of rotation x Q 
for reasons to be made clear in 
the next section. 


Combined angle of attack 
and pitching velocity .- Since 
the pressure coefficient 
(eq. (55)) is proportional to 
the square of the normal com- 
ponent of velocity, there arises 
the possibility that the forces 
and moments due to a motion 


involving both pitching and 
angle of attack may not be 
treated separately, as in the 
case of the usual linear analy- 
sis. Therefore, we begin by 
investigating the forces due to 
a combined motion. Further, for 
the same reason, it is not imme- 
diately evident that the trans- 
fer relations given previously 
(eq. (5^)) remain valid, so that 
rather than choosing a particular 
axis of rotation, we shall con- 
sider the body to be pitching 
about an arbitrary axis, x Q . 
Referred to a fixed system of 
coordinates, the situation is 
as illustrated in sketch (e). 

In the body system of coordi- 
nates, the components of stream 
velocity in the axial, radial, 
and azimuthal directions are. 


u = V - qr sin w 

v = V tan a sin <*> + q(x - Xo)sin u) 
w = V tan a cos u + q(x - XqJcos u 



(57) 
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The components u and v are in turn resolved into the component % 
normal, to the body surface by the relation, 

Vjj = u sin 0 - v cos 0 ( 58 ) 


so that, 

Vjj = V(sin 0 - tan a sin co cos 0) - q. sin w[(x - Xq)cob 0 + r sin 0] 


(59) 

Form i ng the pressure coefficient according to equation (55) > substituting 
into equation ( 5 6), and integrating once, we get. 


C N 


2Cq 

= - | J R(x)£-A(x)cos 1 % + B(x)^~ - i sin 2 o>u + |) - 


— G(x)cos (^(sin 2 ^ + 2) dx 


Cx = f/ R(x)tan 0 |a(x)^o% + ^ - B(x)cos + 

2G(x)(~ - i sin 20 ^ + dx 

r *u r 

Cm. = — / [(x - Xo)R(x) + tan 0 R 2 (x)] -A(x)cos Uu + 

o 

B(x)(^ - i sin 2 (j^ + ^ ^ G(x)cos o> ll (sin 2 a) u + 2)Jdx 


>( 60 ) 


where 
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A(x) 

B(x) 

G(x) 


■ 2 sin 2 0 


= -k sin a[( — -y— ^cos 0+| sin 0J - 2 sin 20 tan a 


■ &) [? 


(x - Xn) 

sin 0 + • — r- ■ — cos 0 | + tan a cos 0 


2 


) 


} 


When the entire body experiences compression flow., the limit of 
integration x u may be replaced by 2, and ^ by 5 >, 'whereupon equa- 
tions (60) simplify to. 


% - ~ sf R ( x ) B ( x ) dx 


C x = / E(x)tan 0 [A(x) + G(x)]dx 




(61) 


Cm = ikf “ Xo)H(^c) + tan 0 R a (x)]B(x)dx 


The latter equations apply so long as the inequality given below is 
satisfied over the entire body: 


t? 

tan” 1 — 


tan a + ~ (x - Xq) 


1 - qR/V 


< 0(x) 


(62) 


For bodies with continually growing cross sections at small angles of 
attack (essentially, a < 0(x)), the inequality (62) will in fact be 
satisfied; attention in this case can, therefore, be confined to the 
simpler set of equations (6l). It is to be noted in equations (6l) that 
both the normal force and pitching moment depend only on B(x) . Then, 
since a must be small by virtue of equation (62), the small-angle 
approximation becomes valid, so that B(x) may be written as the sum 
of two terms, each linear in q. or a alone. Therefore, (% and Cm are 
expressible in terms of stability derivatives in the usual way; that is. 
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Thus, comparing terms with the transfer relations, equations (54), it 
can he seen that under the same conditions for which equation ( 63 ) is 
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valid, the transfer relations axe also valid, for the integrals in the 
last three of equations ( 64 ) correspond to the stability derivatives 


‘S) 

equation (6l). 


'V 


and Cm , as may he verified by letting 
° o 


Xq be zero in 


Stability derivatives for large values of a and q .- For larger 
values of a and qZ/V, for which the inequality ( 62 ) fails, it is nec- 
essary to return to equations (60). In this case, it is still possible 
to write a set of equations analogous to equations (63), although the 
definitions of the stability derivatives must be revised somewhat. Usu- 
ally, in a stability analysis, one is interested in the small deviations 
from some equilibrium condition. This may be represented in equations 
(60) by letting a and q be. 


a = ctrp + Ax 
q. = qrp + Aq 


( 65 ) 


where cxj, and q^, are the angle of attack and pitching velocity corre- 
sponding to the equilibrium (trim) condition, and Ax and Aq are small 
deviations in these quantities. Correspondingly, the normal-force and 
pitching-moment coefficients may be written. 


Cm = Cm(otgi, qrp) + ACm I 




(66) 


where Cjj(otrp, q T ) and C m (a T , q^,) are equations (60) with a and q replaced 
by (Xj, and q^,. Then, since Act and Aq are small, use of the small-angle 
approximation is again permissible. Hence, retaining only first-order 
terms in Ax and Aq, we may write the changes ACjj and ACm as. 


ac n 


\2Ao/ 

Ax- 


Ax + 




*\ 


> 


(67) 
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The quantities in parentheses are the stability derivatives , corresponding 
to the analogous terms in equations ( 63 ). Now, however, they are to be 
evaluated near the equilibrium condition, and their values will depend on 
the angle of attack a.j, and pitching velocity q^ corresponding to that 
condition; that is. 


/So N \ p r* 1 , , 

c N a = = - g J R(x)H(x, ofcj,, qjp)dx 


Aa-»o 


SCw 
^ V 


2 r** 

“ s / R ( X ) J ( X > 


Aq.-^o 






^ = = SlJ [( x " x o)R(x) + tan e H 2 (x)]H(x, c^, q^Jdx 

Aa-^-o o 


> (68) 


c mq = 


dcA 2 n 

— = si I [ ( x - XoiRCx) + tan 9 R 2 (x)]j(x, cup, q. T )dx 


V 


Aq.-> o 


The quantities H and J are, from equations ( 60 ), 


H = (^-i Bto 2 ^ + *)g : .|co 8 Msi n ^ + 2 )^ 


J - (t - I 8in 2Uu + i) 008 ^^ + 2) 


> (69) 


,/AqA 

\v) 


t) 


wherein 
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2(— )[| sin 9 + ( ^p ^cos + tan cx^- sin 28 + g( ^ — ~^ cob 2 Q J 


A few remarks may be in order regarding the use of equations (68). First, 
it can easily be verified that the transfer relations, equations (5*0 , are 
still valid, so that for a given trim condition, it is only necessary to 
compute equations (68) for a single convenient axis position. Also, it 
should be noted that equations (68) ^actually contain (64), as may be veri- 
fied by putting = 0, = ^, and x^ = l in equations (68). 

Second, with regard to the dependence of the stability derivatives on q^ 
and orp: For any practical case it is not conceivable that the pitching 

velocity parameter q^l/v can ever became very large; hence, the depend- 
ence of the stability derivatives on q^ is probably not significant. On 
the other hand, the equilibrium angle of attack cc^ can conceivably be 
very large; in this event, the stability derivatives as evaluated by equa- 
tions (68) can differ significantly from those evaluated according to 
equations (64) and should be used in their place. 


Uniform vertical acceleration .- The impact theory, when applied to a 
uniformly accelerating motion, gives zero for the force and moment pro- 
portional to the acceleration parameter, d l/V. This result is to be 
expected in view of the following considerations. It is known (ref. 10) 
that is closely related to the build-up in normal force that occurs 

following a step change in angle of attack; more precisely, it is propor- 
tional to the area enclosed by the indicial response curve and the steady- 
state ordinate of the indicial curve. Now, implicit in the development 
of impact theory is the assumption that the pressure response to the impact 
of each flow particle is instantaneous. As a consequence, the indicial 
response in normal force to a step change in angle of attack is itself a 
step, whence it follows that the area proportional to Cw. is zero. Like- 

CG 

wise. Cm. is proportional to the area between the indicial pitching-moment 
ct 
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variation and its steady-state ordinate. Again, by impact theory this 
variation is a step,' so that, for the same reason, C^. is zero. 

Stability derivatives for a cone .- F inally , in order to complete 
the set of stability derivatives obtained earlier by the potential 
theories , the stability derivatives for a cone as derived from, impact 
theory are presented below. Calculations are based on equations (64) 
with x Q = 0 , and apply only to the case wherein the independent varia- 
bles (a, d, q.) approach zero. 



2 

1 + T 2 






4 

3 


> 


Cm^ = - (l + t 2 ) 


(70) 


^d = Si = ° J 


It is of interest to note (comparing eqs. ( 70 ) with eqs. (40), (45), 
and ( 53 )) that the centers of pressure for the angle-of -attack and 
pitching cases as derived from impact theory are identical to those 
derived from the potential theories. 


DISCUSSION OF RESULTS 


In order to indicate the nature of the results as obtained from the 
theoretical methods developed herein, numerical calculations have been 
carried out for two cones having semivertex angles of 10° and 20°. 

Results for the variation with Mach number of the stability derivatives 
Cw , Cjt , and Cjj. are shown in figures 1 to 3* The pitching-moment 

M i-v CL 

variations are not shown since in all cases they are simple multiples of 
the normal-force results (see eqs. (4o), (45), ( 50 ), (53), and ( 70 )). The 
curves obtained from potential theory have been terminated at the low end 
of the Mach number scale slightly below the Mach number for which the bow 
wave detaches (ref. ll) , and at the high end, at the Mach number for which 
the Mach cone lies on the body surface. 




! 
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As noted before, it is believed that tbe results for tbe stability 

derivatives due to pitching velocity and to vertical acceleration, 

obtained by use of a first-order and a combination first- and second- 

order potential theory, are comparable in accuracy to the results for 

Cw obtained by the corresponding theory. Therefore, in order to esti- 
oc» 

mate the Mach number range in which the results should apply, consider 
first the results for Cjj a (fig. l) . Here, a comparison is made between- 

the various approximate results and the exact numerical result as obtained 
from reference 2. The curves labeled "first-order, exact Cp" and "first- 
and second-order, exact Cp" are those obtained from calc ulat ions based on 
equations (kO) and (53) • Tbe curve labeled "first-order, approximate Cp" 
is the result, when using the first-order potential solution, of reta ining 
only the first-order term in the expansion of the pressure relation, equa- 
tion ( 28 ) . It is clear from examination of the results that there iB a 
significant improvement in accuracy, even in the first-order solution, if 
the pressure relation is not approximated. Even so, however, only the 
"hybrid" solution can be said to be applicable throughout the Mach number 
range for both cones. Therefore, in the subsequent results, figures 2 
and 3, it is to be assumed that only the hybrid solutions are representa- 
tive of the exact variations for all Mach numbers within the limi ts, the 
Mach number for bow-wave detachment, and the Mach number corresponding 
to 0t = 1. Also shown in figure 1 are the results for Cjf a as obtained 
from the Newtonian impact theory (eqs. ( 70 )). These results are useful as 
a guide for estimating the amount by which the hybrid theory tends to over- 
estimate the magnitude of the normal- force coefficient at the higher Mach 
numbers. A fairing of the hybrid-theory result into the impact-theory 
result is easily accomplished and may serve partially to compensate for 
this tendency. 

Consider next the results for Cjj (fig. 2). It is noted that the 

Q. 

curves for the various approximations are in the same relation to one 
another as were the curves for the corresponding approximations to Cjy^. 

Also, as in the previous case, the hybrid-theory result approaches the 
result obtained from impact theory in a manner which permits a judicious 
fairing of the two. 

Finally, consider the result for Cjj-. (fig. 3) • The interesting 

CL 

point here is that, in contrast to the case of a wing, Cjj^ as obtained 

from potential theory iB positive throughout the Mach number range. This 
fact Implies that at leaBt for axes ahead of the center of loading due 
to d, the damping moment + Cm^ cannot be destabilizing, since the 

destabilizing contribution can only arise from a negative value of C^. 
Actually, it can be shown that Cm. + Cm. is not a destabilizing damping 

3. GG 

moment for any position of the axis by operating with the transfer rela- 
tions (eq., (54)) as was done in reference 10. It is also worthy of note 
in figure 3 that Cjf. vanishes when the Mach cone lies on the body 
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surface. This result is consistent with that obtained from impact theory, 
and, as mentioned previously, implies that at this and greater Mach numbers 
the indicial normal-force response to a step change in angle of attack is 
itself a step. 


Ames Aeronautical laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 21, 195 6 



APPENDIX A 


DERIVATION OF PRESSURE COEFFICIENT 


We derive here the pressure relationships used in the potential- 
theory analysis for each of the three motions considered herein. 

Consider a body ■which moves past an x-l, y x , z x coordinate system 
fixed with respect to still air. Then from Bernoulli’s equation for 
unsteady irrotational flow, one has (ref. 12 ) 




1 /_ 2 . _ 2 . _ 2 \ . 7 P 7 Po 

- (<Py_ + 9y 1 + <P K . ) + i — r - = - — T “ 


7 “ 1 P 7-1 Pr 


(Al) 


Since the flow iB isentropic. 


so that. 


— = constant 
P 7 


p 7 IZ i 

P _ P o p 7 


P Pr 


(A2) 


Substituting equation (A2) in (Al) and solving for p. 


7-3- 7-1 

7 ~T~ 

P = Po 


{ x ■ ^ if K * t (v 2 + v 2 - % £ )] } (A3) 


But ^-2 = a Q 2 , where a Q is the speed of sound in still air. Then 
p o 


- Po i 1 ' K + 5 (% E + V + %")] } 


7_ 

7-3- 


and 


P - P, 


7_ 

7-3- 


2 - {l - [<P tl + | (-PX3. 2 + ^ + ? Zl 2 )] } -1 


(AU) 
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A reference velocity V is now introduced; we define it to be the velocity 
of the body in its axial direction, since this will be a constant for all 
three motions. Then setting 


? ° “ 7M 2 90 (“ 


' *0 


ipV 2 ) 

2 H o J 


we have. 




This equation refers to axes fixed with respect to still air; we next 
consider the transformations to body axes for each of the three body 
motions . 


Sinking With Uniform Vertical Velocity 


As can be seen on sketch (f), 
the fixed and moving axis systems are 
related by, 

X = x x + Vt-L 

y = yx 


z = z x + Vtjtan a 


t = t : 


Then taking derivatives according to. 



Vt, tana 


V tana 


Sketch (f) 


'Px.j^ - Px 

'Pyi = Vy 
Vz x = ^z 

dx dz 

’v = ^ 5ti + 5ti + ** 


equation (A5) becomes, in the body x, y, z coordinate system. 
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7 


_2_JTi_: 

ZjlI tan 


7 M 2 IL 

2 “X v + V 

a+ v 2 + ? )\ x / 


(A7) 


F inally , converting to cylindrical coordinates by tbe 


transformations , 


x = x 

y = r cos a) 
z = r sin gj 


(A8) 


and letting cp^ be zero for steady flow, we have for Cp in tbe moving 
axis system. 



Pitching With Uniform Angular Velocity 



Sketch (g) 
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(All) 


Changing to cylindrical coordinates by equation (A8) , and letting cp-^ be 
zero for steady flow. 



Sinking With Uniform Vertical Acceleration 



The pressure equation (A5) then becomes. 
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Cp = 2 r fi . 2-1 + 3{t + !£!■ + 

^ 7M 2 XL 2 \ V V 2 V v 2 / J 


2L 

7-1 


x } 

(A14) 


Converting to cylindrical coordinates by equation (a8 ) , 


Cp = Nl - 2Ll_i rffe + + 2dtfe cosu+y sin w) 

^3 7M 2 V L 2 L V V 2 \rV 7 / 


(Px 2 + q^ 2 + (q>oj/r) ; 


-»7-x 


: ]J 


- 1 


(A15) 
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